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Abstract 

Considering the most general one-species reaction-diffusion processes 
on a Cayley tree, it has been shown that there exist two integrable models. 
In the first model, the reactions are the various creation processes, i.e. 
oo — >• #0, oo — > •• and o» — >• ••, and in the second model, only the 
diffusion process •o o» exists. For the first model, the probabilities 
Pi{m;t), of finding m particles on l-th shell of Cayley tree, have been 
found exactly, and for the second model, the functions Pi{l;t) have been 
calculated. It has been shown that these are the only integrable models, 
if one restricts himself to L 4- 1-shell probabilities P(mo, mi, • • ■ , m^; t)s. 

1 Introduction 

The integrable reaction-difFusion processes have been investigated by various 
methods on one-dimensional lattice. A class of these models are characterized 
by a master equation with appropriate boundary conditions. These boundary 
terms, vifhich determine the probabilities at the boundary of the space of the pa- 
rameters, are chosen such that the studying of the various reactions becomes pos- 
sible via a simple master equation. The basic quantity in these models is the con- 
ditional probability P(ai, •• ■ ,aN,xi, ■ ■ ■ , xat; • • • , /^at, ■ • • , ?/jv; 0), which 
is the probability of finding particles ai, • • ■ , a^v at time t at sites xi, ■ • • , xat, 
respectively, if at i = we have particles • • • , /Sat at sites yi, - ■ ■ ,yN, respec- 
tively [Di]. 

Another class of models are those which are solvable through the empty 
interval method. In these models, the main quantity, in the most general case, 
is Ek^n{i) which is the probability that n consecutive sites, starting from the site 
k, are empty. Several examples have been investigated by this method [7Ul3j. 

The crucial property which causes the above mentioned cases to be analyt- 
ically solvable, is the dimension of the lattice which particles move on it. One 
dimensional lattice is clearly the simplest case. 

Considering the higher-dimensional lattice in studying the integrable models 
is obviously a great improvement in theoretical physics. One of the important 
example in this area is Cayley tree. A Cayley tree is a connected cycle- free graph 
where each site is connected to z neighbor sites, z is called the coordination 
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Figure 1: Cayley tree with coordination number z — S. 



number. The Cayley tree of coordination number z may be constructed by 
starting from a single central node (called the root of the lattice) at shell ^ = 0, 
which is connected to z neighbors at shell I — 1. Each of the nodes in shell I > 
is connected to (z — 1) nodes in shell I + 1. This construction may stop at shell 
I = L, or continue indefinitely (Fig.l). The interior of an infinite Cayley tree is 
called the Bethe lattice. Due to distinctive topological structure of Cayley tree, 
the reaction-diffusion processes on this graph may be integrable. 

The diffusion-controlled process of cluster growth, introduced by Witten and 
Sander [M] , has been studied on Cayley tree in [15] , and in [16] , the two-particle 
annihilation reaction for immobile reactant has been studied on Bethe lattice. 
Diffusion- limited coalescence, A + A ^ A, and annihilation, A + A |17) . 
and random sequential adsorption [TH] have also been studied on Cayley tree. 
Also the reaction-diffusion processes on Cayley tree, solvable through the empty 
interval method, have been studied in [19j. 

The present paper is devoted to the study of some integrable reaction- 
diffusion processes on Cayley tree. Denoting the number of particles on l-th shell 
by mi, we seek the situations in which the probability P(mo, mi, m2, • • • ,mL;t) 
can be calculated exactly. Zmax = L is the last shell of the Cayley tree. By con- 
sidering the most general reactions on Cayley tree in section 2, it is shown 
that we must restrict ourselves to three distinct creation reactions so that 
a subclass oi L + 1-shell probabilities P{mo,mi ■ ■ ■ ,mL;t), i.e. Pi{m;t) = 
P(0,0, ••• ,0,TOi = m, O,--- ,0;i), becomes solvable. The exact solution of 
Pi{m;t) is obtained in section 3. Finally in conclusion section it is shown that 
the diffusion process can also lead to an integrable model, if one restrict himself 
to one particle probabililities Pi{\; t), and if an extra trapping reaction exists at 
the central node. An ideal spherical trap, surrounded by a swarm of Brownian 
particles, is a fundamental model which had been first presented by von Smolu- 
chowski [22] , and has been generalized to a mobile trap in [21] . Also the energy 
transfer process in dendrimer supermolecule on the Cayley tree with a central 
trap has been discussed in [35]. These two models are only integrable models 
in this context. 

It must be added that the procedure introduced in this paper can be easily 
extended to the situations in which the number of occupied shells is greater than 
one, but there must exist at least one empty shell between any two of them. 
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Figure 2: The source terms of configuration (a) of (0, 1, 2) occupation numbers. 
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Figure 3: Tlie source terms of configuration (b) of (0, 1, 2) occupation numbers. 



2 The integrable models 

The most general reactions of single species models with nearest-neighbor in- 
teractions are 

1 : oo — > #0, ri 4 : #0 — > o», n 

2 : oo — > ••, r2 5 : "o oo, rs 

3 : o» — > ••, ra 6 : •• — ?• oo, rg 

7 : •• "o, rr (1) 

where an empty (occupied) site is denoted by o(»). The r^s are reaction rates 
and there is no distinction between left and right. Our main goal is to find 
the situations which lead to integrable models. By integrability, we mean the 
possibility of the exact calculation of probabilities -P(mo, mi, • • • ;t). 

The necessary condition for achieving this goal is that the occupation num- 
bers of shells (the number of occupied sites of every shell) must be the only 
parameters needed to characterize the configurations. Clearly this is not the 
case if two adjacent shells are occupied. For example consider P(0, l,2;t) and 
the reaction (1) of eq.(l). In Figs. (2) and (3), the source terms of two different 
configurations (a) and (b), with same occupation numbers (0,1,2), are shown. 
As it is seen, the configuration (a) has two source terms, while the configura- 
tion (b) has three. The number behind each configuration is its multiplicity, 
i.e. the number of ways which this configuration can lead to the desired state. 
This shows that if two adjacent shells are occupied, it is not possible to write 
the same evolution equation for different configurations of a given occupation 
numbers. Therefore, we only consider the situations in which there is only one 
non-zero occupation number, i.e. — mdu. As mentioned earlier, it is possible 
to extend our procedure to the situations in which there is at least one empty 
shell between any two occupied shells. 

Considering the probabilities Pi{m; t) — P(0, 0, • • • ,0,mi — m, 0, • • • ,0;t), 
we must yet check whether all the reactions of eq.([T]) are acceptable or not. By 
the term acceptable, we mean to have two properties. First, the determination 
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of occupation number m is sufficient for liaving a unique evolution equation, 
irrespective of particles' distribution. Second, all the terms of evolution equation 
are expressible in terms of Pk{n; t)s. Now it can be seen that this is not the case 
for reactions (4)- (7) of eq.((T]). This is because in all these cases, the source terms 
of P/(m; t)s are 2-shell probabilities P{mi-i,m,i; t) and P(rni, mz+i; t), which can 
not be expressed in terms of Pfc(n; t)s. For these reactions, the master equations 
are: 

d 

Pi{m; t) = aiP{mi^i = 1, m; = m — 1; t) + biP{mi — m — 1, m/+i = 1; t) 



dt 



~SiPi{m-t) (for r^n), (2) 

d 

—Pi{m;t) = a5P{mi^i = l,mi = m;t) + b5P{mi = m,mi+i = l;t) 

+C5Pi{m+l;t)^S5Pi{m;t) ( for r = rs), (3) 

d 

Pi{m;t) = a(iP{mi^i = l,nii = m + l;t) 



dt 



+beP{mi = m + l,mi+i ^ l;t) (for r = re), (4) 



and 



d 

—Pi{m;t) = a7P{mi^i = l,mi = m;t) 



dt 



-b-jP{mi^m,mi+i^l;t) (for r = rr), (5) 



respectively. The parameters ai,bi,Ci and Si are some constants. The case 
m = 1 is an exception which will be discussed in section 4. 

So the integrable model, through the Pi{m;t) probabilities, is a model in 
which the allowed reactions are reactions (l)-(3) of eq.(IT]). The evolution equa- 
tion of P/(m; t) on a Cayley tree with coordination number z and L + 1 shells 
then becomes (for ^ = 0, 1, 2, • • • ,L~1) 

d 

— Pi{m;t) = {ni - m + l)zriPi{m - l;t) 



dt 



In above equation 



[2(y~^ nil — mz)ri + (^^ ni — mz)r2 + mzr^]Pi{m\ t).(6) 



, z(z- for l> 1 

"'-^ 1 for / = 



is the number of sites of the ^-th shell. The first term in the right-hand site of 
eq.® is the source term of reaction (1), and 2(^p^-^ n// — mz)riPi{m; t) is its 
sink term. This is because X]i'=i ~ "^•^ the number of empty pairs (o o) 
of a Cayley tree with m particles on shell /, which each of them can change to 
either • o or o • with rate ri. The reaction (2) has no source term for P/(m; i), 
and its sink term is clearly (X]i'=i ~ mz)r2Pi{m;t). Finally the reaction 3 
does not have source term, and its sink term is mzr^Pilm] t) because each of the 
neighbors of each particles of l-th shell can be changed via o» — >• •• reaction. 

For I — L, since there is no upper shell, the number of neighboring sites of 
each occupied site is one, instead of z. Therefore for PL{m]t), the evolution 
equation is same as eq.®, when z is replaced by 1. 
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3 Solutions of the master equation 

In this section we try to solve the master equation 
d 

—Pi{m\t) = {ni-m + l)zriPi{m-l]t) 

-~[{r:i-r2-2ri)mz + {2ri+r2)N]Pi{m-t) {I < L), (8) 

with 

L 

iV^^n^s (9) 

i'=i 

in two cases — r2 — 2ri ^ and — r2 — 2r\ — 0. For I — L, we must set 
z = 1 in eq.®. 

3.1 The case ra 7^ r2 + 2ri 

To solve eq.®, we use a recursive method. We first consider the case m = 0, 
where eq.® is reduced to 

^P,(0;t) = -iV(2ri+r2)P,(0;t), (10) 

with solution 

P/(0;i) = P,(0)e-^(2n+r-2)t^ ^^^^ 

Pi(0) = Pi{Q]t = 0) is the probability of finding no particle in shell I &tt = {). 
Inserting eq.((TT]) in m = 1 case of eq.®, results in Pi{l]t) as follows 



P.(l)- -^^i^KO) 



g-[z(r3-r2-2ri)+Ar(2ri+r2)]t 



niriPijQ) ^„jv(2n+r,)t ^3) 



Continuing this method, one can deduce the following general expression for 
Pi{m-t): 

m k f -w / \ rn-j 

P(rn-t) = W ^ (-l)fe-^- 

^^^(n,-m)!(m-A:)!(fc-j)! ^r3-r2-2rj ^ ^ 

xPi(j)e-[^'^(''^-'^^-2ri)+Ar(2ri+r2)]t (/<i)^ (I3) 

in which, as we will show, Pi{j) = Pi{j; t = 0). 

It can be seen that for m — and m = 1, eq. p^ leads to eqs.((TT]) and (HU, 
respectively. To prove that ea. dT^ is the solution of eq.®, we insert it in eq.®. 
The summations in the left-hand side and in the second term of the right-hand 
side of eq.® are both Y,T=o I^j=o- '^'^ write them as YlT=o + YlT=o Z]j=0' 
then it can be seen that the X^JLo terms in both sides are the same, and therefore 
they are cancelled. So in all three terms of eq.®, the summations have the same 
form X^I-Lo^ X]j=0' 'which after some simple calculations, it can be shown that 
the remaining terms are cancelled, which proves the expression (1131) satisfies 
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the master equation ([S]). Note that because of the (ti; — m)! factor m the 
denommator of eq. p^ . Pi{m;t) satisfies 

Pi{m>ni;t)=0, (14) 

which shows the correct behavior of eq. dT51) . 

It is also necessary to prove the physical interpretation of Pi{j) as Pi{j; t = 
0). To see this, we consider ea. (fT3)) at t = 0: 

mm/ s m—j , 

■^^^ V'"3 - f2 - 2riy (r7,;-m)! 

^(m-fc)!(fc-j)!' ^^^^ 
in which we use the equality: 

m k mm 

5]^A(j,fc)=^^^0-,fc). (16) 

k=0 j=0 j=0 k=j 

Using k' = k — j, eg. dTSl) becomes 

m , \ rn-j , .s. 

^ V'^s - ''2 - 2ri/ (nj-m)! 

^ k'Um~j~k'y: ^ ' 

k'=0 ^ ■' ' 

From binomial expansion 

k'=0 ^ ' 

one finds for a = 5 

= n!a"V ,,;. (19) 
^ fc' (n - fc') 

So the last summation of ea. (fT7|) is zero for all m — j 7^ 0, or j = 0, 1, ■ • • , to — 1. 
The only remaining term is j = to, which results in 

Pi{m-t = ^) = Pi{-m). (20) 



This completes the proof of eq. p3|) as the exact solution of the master equation 
with Pi{j)s as the initial values of probabilities. The probability PL{m\t) 
is found from (fT3)) . by taking I — L and z = 1. 

Eq. p^ gives the probability of finding to particles at time t on shell /, when 
all other shells are empty, if we begin with any number of particles on shell I at 
t = 0. Of course for a specific initial condition, i.e. Pi{j',t = 0) = Sjjg, only one 
term of summation X]j=o survives. 
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3.2 The case = r2 + 2ri 

For the case ^3 = r2 + 2ri , the master equation ([5]) becomes 
d 

—Pi{m-t) = {ni-m+l)zriPi{m-l-t)^NnPi{m]t) il < L). (21) 

Using the method of the previous case, the solution of above equation is found 
to be 

™ ( — '\\ 

P^i^^.t) = e-''-^'Y.l \u ,, Pi{j){riztr-' (KL). (22) 

•■f— ^ [m ~ j)\[ni — my. 



It can be easily shown that the solution ([22)1 satisfies (1211) and has the desired 
properties Piijn > ni;t) = and Pi{m;t = 0) = Pi{m). For I = L, the solution 
is found from ea. ip^ . by taking I — L and z = 1. 

4 Conclusion 

To study the integrable reaction-diffusion processes on a more-than-one-dimensional 
lattices, the most general interactions ([T]) have been considered on a Cayley tree 
with coordination number z. It has been shown that among the probability 
functions, the probabilities Pi{m;t)s are independent of distribution of parti- 
cles on various shells, and may lead to integrable models if we consider two 
situations. The first one is the creation-reactions (l)-(3) of eq.(IT|), with master 
equation ([SJ. The Pi{m;t)s for two cases ^ r2 + 2ri and — r2 + 2ri 
have been considered in section 3, with the final exact results (IT5|) and (1^^ . 
respectively. 

As is clear from eqs.(l2])-([5]), in the case m — 1, the master equations ©-([S]) 
still contain the 2-shell probabilities, but for diffusion process r = r^, only the 
one-shell probabilities (i.e. one point functions) appear in eq.©, which may 
lead to an integrable model. 

The master equation oi Pi{l;t), with I > 1, for diffusion process on a Cayley 
tree with coordination number z is: 

§-^m;t) = {z- l)Pi-iil;t) + Pi+,il;t) - zP/(l;i), (23) 

in which the diffusion rate has been absorbed in the rescaling of time t. The 
evolution equations of Po{l;t) and Pi{l;t) are 

^^Po{l;t) ^ Piil;t) ~ zPoil;t), (24) 

and 

^Pi{l;t) = zPQ{l;t) + P2{l;t) ^ zPi{l;t), (25) 
ot 



respectively. These equations take a form similar to ea. ([23|) . provided one defines 

F_i(l;0:=0, (26) 
Po(l;i):=0. (27) 
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The first boundary condition is trivially satisfied, since Pi (m; t)s have been 
defined for I > 0. But the second one indicates a trapped reaction at the origin. 
So ea. (P5)) . with arbitrary I, defines an integrable model on a Cay ley tree if 
the particles have diffusion process, and if there exists a trap at the root of 
the lattice. When the reactions are coalescence A + A ^ A and annihilation 
A + A 0, the density of particles in shell i.e. pi{t), has been calculated 
in [T7] for the same situation, that is a trap at the root of a Cayley tree. 

It can be shown that the solution of master equation (j23p with boundary 
condition (P7)) is 

Pi{l;t) - (z - (2(z - ly/H) - /w„(2(z - 1)1/2^)] e-*, (28) 

in which /„ is the n-th order modified Bessel function. 

Finally it must be added that it remains one more case which may lead to an 
integrable model. Looking at eqs.©-©, it is clear that if one takes m = 0, only 
the master equation ([3]) leads to an equation which only consists of one-point 
functions, i.e. 

^Pi{0;t) = a5Pi-i{l;t) + b5Pi+i{l;t) + c5Piil;t) - s5Pi{0;t). (29) 
ot 

It can be shown that S5 = 0. But the point is that the determination of Pi{0; t)s 
depends on the evaluation of Pi{l;t)s, which can not be calculated. This is 
because the master equation of Pi{l;t)s contains the two-point probabilities 
P(l, l;t) and therefore is not closed. Therefore the two considered models are 
the only cases which can be exactly solved. 
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